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ABSTRACT 

Context. Recent studies have started to cast doubt on the assumption that most stars are formed in clusters. 
Observational studies of field stars and star cluster systems in nearby galaxies can lead to better constraints on the 
fraction of stars forming in clusters. Ultimately this may lead to a better understanding of star formation in galaxies, 
and galaxy evolution in general. 

Aims. We aim to constrain the amount of star formation happening in long-lived clusters for four galaxies through the 
homogeneous, simultaneous study of field stars and star clusters. 

Methods. Using HST/ACS and HST/WFPC2 images of the galaxies NGC 45, NGC 1313, NGC 5236, and NGC 7793, 
we estimate star formation histories by means of the synthetic CMD method. Masses and ages of star clusters are 
estimated using simple stellar population model fitting. Comparing observed and modeled luminosity functions, we 
estimate cluster formation rates. By randomly sampling the stellar initial mass function (SIMF), we construct artificial 
star clusters and quantify how stochastic effects influence cluster detection, integrated colors, and age estimates. 
Results. Star formation rates appear to be constant over the past 10'^ - 10^ years within the flelds covered by our 
observations. The number of clusters identifled per galaxy varies, with a few detected massive clusters (M > 10^ M©) 
and a few older than 1 Gyr. Among our sample of galaxies, NGC 5236 and NGC 1313 show high star and cluster 
formation rates, while NGC 7793 and NGC 45 show lower values. We flnd that stochastic sampling of the SIMF has a 
strong impact on the estimation of ages, colors, and completeness for clusters with masses < 10^ - 10^ Mo, while the 
effect is less pronounced for high masses. Stochasticity also makes size measurements highly uncertain at young ages 
(r < 10^ yr), making it difficult to distinguish between clusters and stars based on sizes. 

Conclusions. The ratio of star formation happening in clusters (F) compared to the global star formation appears to 
vary for different galaxies. We find similar values to previous studies (F ^ 2%-10%), but we find no obvious relation 
between F and the star formation rate density (Ssfr) within the range probed here (Ssfr 10~^ — 10~^ Mq yr~^ 
kpc~^). The F values do, however, appear to correlate with the specific U-band luminosity (Ti:(U), the fraction of total 
light coming from clusters compared to the total U-band light of the galaxy). 

Key words, galaxies: Individual - NGC 5236 galaxies: Individual - NGC 7793 galaxies: Individual - NGC 1313 galaxies: 
Individual - NGC 45 galaxies: Star clusters - galaxies: Star formation - galaxies: Photometry 



1. Introduction 

It is often assumed that stars can be formed either in the 
field of a galaxy as single stars or in a group of stars (cluster) 
that formed from the same molecular cloud at the same 
time. However this v iew has been recently questioned by 
iBressert et alJ (|2010l ), who challenged the idea that field 
and cluster formation are actually distinct modes of star 
formation. 

O wing to dyna mical and stellar evolution clusters dis- 
rupt ([Spitzerl [19871 ) and the stars become members of the 
field stellar population. While it is commonly assumed that 
most (if not all) stars formed in clusters (e.g. iLada fc Ladal 
l2003l : [Porras et al.i .2003, in the solar neighborhood), the 
amount of star formation happening in those clusters that 
remain bound beyond the embedded phase is still uncer- 
tain. In any case, if clusters are to be used as tracers of 



* Table 10 is only available in electronic form at the CDS 
via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via 
|http: / / cdsweb.u-strasbg.fr / cgi-bin / qcat? J / A+ A / , 



galactic star formation histories, it is of key importance to 
understand what fraction of star formation is happening in 
long-lived clusters and whether or not this fractio n corre- 
l ates w ith other host galaxy parameters. Following iBastianI 
(I2008D . we refer to this fraction as F. 

Estimating F is not straightforward. Apart from dif- 
ferences arising at the time of formation, cluster disrup- 
tion will also affec t the detect e d nu mber of clusters of 
a given age (r). ILada fc Ladal (|2003t ) estimate that be- 
tween 70% to 90% of the stars in the solar neighborhood 
form in embedded star clusters, while only 4-7% of these 
clusters survive for more than about 100 Myr. Similarly, 
[Lamers & Gieles (2008) estimate an "infant mortality" rate 
of 50% to 95%, based on a comparison of the surface den- 
sity of open clusters and the star formation rate near the 
Sun. By studying t he UV flux in and ou t of clusters in the 
galaxy NGC 1313, "Pellerin et al.' (2007) suggest that over 
75% (between 75% and 90%) of the flux is produced by 
stars in the fleld, concluding that the large number of B- 
type stars in the fleld of the galaxy could be a consequence 
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of the (high) infa nt mortality of cl u sters. For the Small 
Magellanic Cloud, iGieles fc BastianI (|2QQ8[ ) estimate that 
optically visible, bou nd clus t ers ac count for 2%-4% of the 
star formation, while "Gieles| (|2Q1Q[ ) estimates this fraction 
to be in the range 5%-18% for the spiral galaxies M74, M51, 
and MlOl. However, most of these studies could not distin- 
guish between scenarios in which a large fraction of stars 
initially form in clusters that rapidly dissolve or whether 
there is a genuine ^^fie ld" mode o f star formation. Studying 
the Antennae galaxy, iFalll (|2QQ4l ) estimates that 20% (and 
possibly all) stars were formed in cluste rs. Using a larger 
sample of galaxies, Goddard et al.' ("2010) find a power-law 
relation (F cx Sgpj^) between the fraction of stars forming 
in clusters that survive long enough to be optically visible 
and the star formation rate density of the galaxy (Ssfr)- 
Their data set covers different types of galaxies, from ir- 
regulars (i.e. LMC, SMC, and NGC 1569) to grand design 
spirals (i.e. NGC 5236). The Ssfr of these galaxies vary 
from 7 X 10"^ to - 700 x 10"^ Moyr^^Kpc-^ Their F vs. 
^SFR relation, however, is based on somewhat heteroge- 
nous data with different mass- and age ranges, which do 
no t come from the same observations (see details in Sect. 4 
of iGoddard et alll2010f ), although the authors do attempt 
to homogenize the sample by normalizing the cluster sam- 
ples to a common mass limit. 

The actual definition of the phase called infant mortality 
is somewhat ambiguous in the literature. Early disruption 
due to rapid gas expulsion may only take a few Myrs, but 
the term has also been used to describe mass-independent 
disruption, meaning that clusters lifetime is independent 
of mass over a much longer time span. In the latter case, 
the "infant mortality rate" (IMR) refers to the fraction of 
clusters that are disrupted per decade of age. We prefer to 
simply use the term "mass-independent" disruption (MID) 
in this case. For MID, the IMR is related to the slope a of 
the age distribution, dN/dr oc of a mass limited clus - 
ter sample as a = log:(l - IMR) (IWhitmore et al.l l2QQ7t ). 
For example, 'de G riis fc GoodwinI (|2008^ found that for 
the SMC the IMR is close to 30% (between 3-160 Myr), 
while the logarithmic age distribution of clusters in the 
Antennae galaxies is about fi at (a ~ —1), indicating an 
IMR close to 90% (|Fallll2004D . assuming that the star for- 
mation rate has been about constant over the past few 10^ 
years. On theoretical grounds, the time scale for the (grad- 
ual) cluster disruption is expected to be mass-dependent, 
owing to tidal sh ocks and eva poration that follows early 
gas expulsion (e.g. lGieles et al.| [2006), assuming there is no 
strong relation between cluster mass and radius. In this 
description, the dissolution time t^is of a cluster scales 
with cluster mass as t^is = t^j M/ 10^ Mc:))^ , where t^ is 
the lifetime of a 10^ M^;^ c luster (see iBoutloukos fc LamersI 
l2003l : iLamers et al.l |2005[ ) . The time scale on which clus- 
ters dissolve may also depend on external factors, such 
as the tidal field strength, density of molecular gas, pas- 
sages near /through gia nt molecular clouds, or t hrough spi- 
ral arms, etc. (see e.g. iGieles et al1l2006L 120071 ). This sce- 
nario attempts to compile in one single formula all the pos- 
sible p rocesses that affect cluster disruption. See ILamersI 
(|2009f ) for a description of the different models for cluster 
dissolution. 

Determining the extent to which cluster dissolution is 
a mass-dependent process has turned out to be difficult. 
Estimations of cluster parameters based on observations 
are affected by stochastic effects, degeneracies, and obser- 



Fig. 1. Galaxies studied in this paper. Top left: NGC 5236; 
top right: NGC 7793; bottom left: NGC 1313; and bottom 
right: NGC 45. Red lines represent the pointings covered 
by the HST/ACS, while the blue lines represent the point- 
ings of the HST/WFPC2. Images were taken from the DSS 
archive using Aladin software. 



vational uncertainties. For example, iMaiz ApellanizI (|2009f ) 
used Monte Carlo simulations to estimate how stochastic 
effects coming from the random sampling of the stellar 
initial mass function influence the determination of ages 
and masses, which a re de rived from broadband photome- 
try. iPiskunov et al.l (|2009l) show how the consideration of 
the discreteness of the stellar initial mass function (IMF) 
can explain features observed in the color-age relation and 
can improve the flt between models and observations. They 
conclude that the large number of red outliers can be ex- 
plained as a systematic offset coming from the difference be- 
tween discrete- and continuous-IMF at low masses (Mc=10^ 
M0) and young ages {\og{r)[yr] ~7), reaching up to ~0.5 
magnitudes, and decreases down to '^0.04 magnitudes at 
higher masses (Mc=10^ M©). 

To estimate fleld star formation histories, a differ- 
ent approach is needed than for clusters, because ages 
canno t in genera l be d etermined directly for individual 
stars. iTosi et al] (|1991[ ) presented a method that takes 
incompleteness, resolution, depth, and observational er- 
rors (among other parameters) into account to construct 
a synthetic color- magnitude diagram (CMD), which can 
be used to estimate the star formation history by com- 
parison with observations . This method has been de- 
veloped further by other authors in the past years, e.g. 
"DolphinI (|1997[ ) and iHarris fc Zaritskvl (|2QQll) . and has 
been used for a lar ge nu mber of galax ies, e.^. SMC. LMC 
(|Harris fc Zaritskvl l2004l l2009D . M31 (|Brown et al.l 120081 ). 
NGC 1313 (ILarsen et al.l[200l . In this series of papers, we 
make use of this method to estimate the fleld star forma- 
tion rates of our target galaxies, which we then compare 
with c luster formation rate s to es timate F. 

In ISilva- Villa fc LarsenI (|2010l hereafter Paper I), we 
presented the tools needed to study and constrain the F 
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value of our set of galaxies, and used NGC 4395 as a testbed 
galaxy. As the second paper in a series, this paper aims to 
estimate T in differe nt environments and comp are it with 
previous work (e.g. iGieiesll2QlQl : ICoddard etldll2QlQt ). us- 
ing the complete set of galaxies. To this end, we took ad- 
vantage of the superb spatial resolution of the Hubble Space 
Telescope (HST) and used images of the galaxies NGC 5236, 
NGC 7793, NGC 1313, and NGC 45, which are nearby, face- 
on spiral galaxies that differ in their current star formation 
rates and morphology. These galaxies are near enough 4 
Mpc) to allow us to disentangle the cluster system from the 
field stars, making it possible to estimate cluster and star 
formation histories separately and simultaneously from the 
same data. 

The paper is structured as follows. In Sect. 2, we present 
a short overview of previous work on our target galaxies, 
related to the present study. The basic reduction and char- 
acteristics of the observations are described in Sect. 3. In 
Sect. 4 we present the photometry procedures applied to 
the data and describe how completeness tests were carried 
out. We also discuss the effect of stochastic sampling of the 
stellar IMF on integrated cluster properties. In Sect. 5 we 
present the results of the estimation of ages and masses of 
clusters, as well as the field star formation histories. We 
also estimate the cluster formation rates and use these to 
determine T values. In Sect. 6 we discuss our results and 
finally, we summarize and conclude our work in Sect. 7. 

2. Dataset overview 

In this paper we describe results for the remaining four 
galaxies in our HST/ACS sample: NGC 5236, NGC 7793, 
NGC 1313, and NGC 45. These four galaxies share the 
properties of being face-on, nearby spirals; however, they 
differ in their morphology, star, and cluster formation his- 
tories. We present the basic properties of each galaxy in 
Ta ble m 

iLarsen fc Richtlerl (|1999) studied cluster populations in 
a set of 21 galaxies, including the four included here. Using 
ground-based multiband {UBVRI and Ha) observations 
they estimated the total number of young massive clusters 
in each galaxy, us ing a magnitude l i mit o f My < —8.5. 
In a further work, Larsen fc Richtlerl (|2QQQ[ ) estimated the 
star formation rate density (Ssfr) and the specific [/-band 
luminosity, Tl(U) = 100 xL(clusters,U)/L (galaxy, U), for 
each galaxy. The Tl{U) was found to correlate with Ssfr- 
Taking Tl(U) as a proxy for the cluster formation efficiency, 
these data thus suggested an increase in the cluster forma- 
tion efficiency with Ssfr- It is worth noting here that the 
SsFR values were derived by normalizing the total star for- 
mation rates, obtained from IRAS far-infrared fluxes, to 
the optical galaxy diameters obtained from the RC3 cata- 
log. Therefore, while these numbers were useful for studying 
trends and correlations, they should not be taken as reliable 
absolute values. 

More recent e stimates of Ssfr have been made 
by IChandar et aD (2010b) for the galaxies NGC 5236 
and NGC 7793, where they found s imilar values to 
ILarsen fc Richtlerl (|2000l ). iHarris et al] (|200ll ) present a 
photometric observation of clusters in the center (inner 300 
pc.) of NGC 5236. Harris et al. find a large number of young 
and massive clusters, consistent with a burst of star forma- 
tion that began around 10 Myr ago, but note that the ap- 
parent absence of older clusters might also be due to rapid 



disruption. IChandar et al.l (|2010cl ) used the new Wide Field 
Camera 3 (WFC3) on HST to analyze the cluster system of 
NGC 5236. They find that luminosity functions and age dis- 
tributions are consistent with previ ous work on galaxies of 
differe nt mo rphological types (e.g. lFallll2003 ). iMora et al.l 
(2003, l2009f ) studied the cluster system for the same set 
of galaxies used in this work, based on the same HST im- 
ages. They present detailed estimates of the sizes, ages, and 
masses for the clusters detected. Mora et al. conclude that 
the age distributions are consistent with a ~80% MID per 
decade in age up to 1 Gyr, but could not make a distinction 
between different models (MDD vs. MID) of cluster disrup- 
tion. In the galaxy NGC 45 they found a large number of 
old globular clusters, of which 8 were spectroscopically c on- 
firmed to be ancient and metal-poor ora et al.ll2008D . 



3. Observation and data reduction 

The five galaxies studied in this series of papers were 
selected for detailed observations with the Advanced 
Camera for Surveys (ACS) and Wide Field Planetary 
Camera 2 (WPFC2) onboar d HST from the work of 
ILarsen fc Richtl^ ([T9991 12OOOI) . The two instruments have 
a resolution of 0'.'05 and 0'.'046, O'.'l for ACS and WFPC2 
(PC,WFs), respectively. At the distance of our galaxies 
(~ 4 Mpc) the ACS pixel scale corresponds to ~ 1 pc. 

Besides NGC 1313, which has three different fields ob- 
served, the rest of the galaxies were covered using two 
pointings (see Fig. [T]). The bands used for the observa- 
tions were F336W(- U), F435W(- B), F555W(- V), and 
F814W(~ /), with the exposure times listed in Table |2l 
The standard STScI pipeline was used for the initial data 
processin g. ACS images were dri zzled using the multidriz- 
zle task (jKoekemoer et al.l 12002") in the STSDAS package 
in IRAF using the default parameters, but disabling the 
automatic sky subtraction. WFPC2 images were combined 
and corrected for cosmic rays using the crrej task using the 
default parameters. 

Object detection for field stars and star clusters was 
performed on an average B, V, and I image, using 
daofind in IRAF for t he stars and SExtractor V2.5.0 
(iBertin fc Arnoutsll 19961 ) for the clusters. Coordinate trans- 
formations between ACS and WFPC2 used IRAF. For de- 
tails we refer to Paper 1. 



4. Photometry 

We review here the procedures for carrying out photometry 
on our data, however, for details, we refer to Paper I. Due to 
the crowding we performed PSF photometry for field stars, 
while we used aperture photometry for the star clusters. 

4.1. Field stars 

With a set of bona-fide stars visually selected in our im- 
ages, measuring their FWHM with imexam, we constructed 
our point-spread function (PSF) using the PSF task in 
DAOPHOT. This procedure was followed in the same man- 
ner for each band (i.e., B,V, and I). The PSF stars were se- 
lected individually in each band, in order to appear bright 
and isolated. PSF photometry was done with DAOPHOT 
in IRAF. 
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Galaxy 


Type^ 


{m-My 




Z 


12 + log{0/H) 


NGC 5236 


SAB(s)c 


27.84 


0.29 


0.008,0.019' 


8.2-8.6' 


NGC 7793 


SA(s)d 


27.6 


0.08 


0.008,0.019^ 


8.57^ 


NGC 1313 


SB(s)d 


28.2 


0.47 


0.004,0.008^ 


8.33^ 


NGC 45 


SA(s)dm 


28.42 


0.09 


0.004,0.008^ 





NASA/IPAC Extragalactic Database (NED}; 
iBresolin etHI ([2009i) : ^ iChandar et al.1 { 
IZaritskv et al.1 (|1994l ) at r 



v^-^^' ^ I ^Q^^ (120091) a n d references therein ; " iSchlegel erall (119981 ): ' 

iChandar et al.1 (|2010bl l:^ I Walsh k Rovl (|1997l ): H^rsen et al.1 (pQQTl ): ^ iMora et al.1 (|2007l ): and 
3 Kpc. 



Table 2. Journal of the observations. 



Galaxy 


Number 


F336W(U) 


F435W(B) 


F555W(V) 


F814W(I) 


RA. 


DEC. 


Date 




of field 


sec. 


sec. 


sec. 


S 


(J2000) 


(J2000) 




NGC 5236 


1 


2400 


680 


680 


430 


13:37:00 


-29:49:38 


2004.07.28 




2 


2400 


680 


680 


430 


13:37:06 


-29:55:28 


2004.08.07 


NGC 7793 


1 


2400 


680 


680 


430 


23:57:41 


-32:35:20 


2003.12.10 




2 


2400 


680 


680 


430 


23:58:04 


-32:36:10 


2003.12.10 


NGC 1313 


1 


2800 


680 


680 


676 


03:18:04 


-66:28:23 


2004.07.17 




2 


2800 


680 


680 


676 


03:18:17 


-66:31:50 


2004.12.18 




3 


2800 


680 


680 


676 


03:17:43 


-66:30:40 


2004.05.27 


NGC 45 


1 


2400 


680 


680 


430 


00:14:14 


-23:12:29 


2004.07.05 




2 


2400 


680 


680 


430 


00:14:00 


-23:10:04 


2004.06.01 



HST zeropoint^ were applied to the PSF magnitudes 
after applying aperture corrections (see Sect. 4.3). The zero- 
points used in this work are ZPb = 25.767, ZPy = 25.727 
and ZPj = 25.520 magnitudes. Typical errors of our pho- 
tometry do not change dramatically from the ones in Paper 
I (see its Fig. 2). 

Having magnitudes for our field stars, Hess diagrams 
were constructed and are depicted in Fig. [2] (each panel 
presents all fields combined for each galaxy). The total 
number of stars varies among the galaxies, all having some 
tens of thousands. Various phases of stellar evolution can 
be identified in the Hess diagrams: 

— main sequence and possible blue He-core burning stars 
at V - / - and -2 < V < -8; 

— red He core burning stars at 1.2 < V — I < 2.5 and 
-2.5 <V < -6.5; 

— RGB /AGE stars, near the detection limit at 1 < y— / < 
3 and -0.5 <V < -2.5. 

The same features were observed for NGC 4395 in Paper 

L 

Overplotted in Fig. [2] are the 50% completeness lines 
(see Sect. 4.5 for details of the completeness analysis). Also, 
red lines enclose the fitted areas that will be used in Sect. 
5 to estimate the star formation histories of the galaxies. 
These areas were selected to cover regions that were clearly 
over the 50% completeness and represent stars in different 
stages of evolution (e.g. main sequence, red He core burn- 
ing). 

4.2. Star clusters 

To detect the cluster candidates we used SExtractor with a 
detection criterion of six connected pixels and a threshold 
of 10 sigma above the background level. The total numbers 



NGC 5236 




-0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 





-0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 -0.5 0.0 0.5 1.0 1.5 2.0 2.5 3.0 

V-l V-l 

Fig. 2. Hess diagram for the field stars of the observed 
galaxies. The dashed white line represents the 50% com- 
pleteness curve. Red lines enclose the fitted areas used to 
estimate the SFH (see Sect. 5). 



www.stsci.edu/hst / acs / analisys / zeropints / #tablestart 



of objects detected in each galaxy are listed in the second 
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column of Table [3l For these objects aperture photometry 
was performed using an aperture radius of six pixels on our 
ACS pointings, corresponding to about two half-light radii 
for a typical star cluster. We used a sky annulus with five 
pixels width and an inner radius of eight pixels. For the 
WFPC2 images, the apertures used cover the same area. 
There is a possibility of having close-neighbor objects that 
contaminate the photometry, whether inside either of the 
aperture radii or the sky annulus. Sizes were n ieasured usin g 
the ISHAPE task in the BAOLAB package (|Larsenlll999l ). 

As mentioned in Paper I, three criteria were used to 
produce catalogs of cluster candidates from the initial 
SExtractor output. 

1. Size: Candidates must satisfy FVKi^MsExtractor > 2.7 
pixels and FI^i^Mishape > 0.7 pixels. These are rather 
conservative size cuts that may eliminate some of the 
most compact clusters, but reduce the risk of contami- 
nation from other sources. 

2. Color. Candidates must satisfy V — I < 1.5. 

3. Magnitude: Candidates must be brighter than my = 
23 {My brighter than —4.6 to —5.4, depending on the 
galaxy distance). 

Since the WFPC2 fields only cover about half the area 
of the ACS fields, some objects will only have three-band 
photometry {BVI)^ while others will have all four colors. 
Objects that satisfy the three criteria listed above are con- 
sidered as star cluster candidates in the rest of the paper. 
However, as found in many previous studies, there is no 
unique combination of objective criteria that can lead to 
a successful detection of bona-fide clusters and no false 
detections. Our cluster candidates were therefore visually 
inspected to determine whether they resemble star clus- 
ters. Based on this, we classified the cluster candidates 
into three categories: Accepted, Suspected, and Rejected. 
Figure [3] presents some examples of each category. In this 
figure, the first row presents the Accepted objects, which 
are clearly extended objects with normal measured sizes 
and magnitudes. The second row presents the Suspected 
objects, where the size/magnitude measurements may be 
affected by crowding, where the shape appears irregular, or 
where the contrast against the background is not strong. 
The last (third) row presents examples of the Rejected ob- 
jects. 

Table [3] summarizes the total number of objects de- 
tected that have size measurements (2^^ column), the total 
number of objects with three-band photometry and have 
sizes over the limits imposed (3^^ column), the total num- 
ber of accepted objects with three- and four band photom- 
etry (4*^ and 5*^ columns), the total number of suspected 
objects with three- and four band photometry (6^^ and 7^^ 
columns), and the total number of rejected objects (8^^ col- 
umn). Shaded areas are the total numbers per galaxy. 

Figure m shows two-color diagrams for accepted plus sus- 
pected clusters with four band photometry (all the fields 
combined per galaxy), corrected for foreground extinction 
with the values presented in Table [H Overplotted is a the- 
oretical track that a cluste r will follow between 4 Myr and 
1 Gy r using Galev models ([Anders fc Fritze-v. AlvenslebenI 
[lop), assuming LMC metallicity and no extinction. We see 
that the clusters generally tend to align with the model se- 
quence, but with significant scatter around it. Below we 
investigate to what extent this scatter may come from 




Fig. 3. Examples of objects that were accepted (first row), 
suspected (second row), and rejected (third row) after vi- 
sual inspection. Stamps are from the first field observed in 
NGC 7793 and have sizes of 100x100 pixels. 



stochastic color variations due to random sampling of the 
stellar IMF. 

4.3. Completeness 

Completeness analysis was carried out separately for field 
stars and star clusters to account for both populations. 

4.3.1. Field stars 

As in Paper I, we created artificial stars using the PSF 
obtained in Sect. 4.1. In the magnitude range between 20 
to 28, every 0.5 magnitudes, we generated 5 images and 
passed each one through the photometry procedures, us- 
ing the exact same parameters as are used for the original 
photometry. A total of 528 stars were added to each image, 
with a separation of 100 pixels (we did not take subpixels 
shifts into ac count). The i mages were created using mksynth 
in BAOLab (|Larsenlll999[ ) and added to the science images 
using imarith in IRAF. To quantify the dependency of com- 
pleteness functions on color, we made use of the near 1:1 
relation between the B — V and V — I colors of stars (see 
Paper I for details). 

Based on our analysis, we found 50% completeness lim- 
its for each galaxy and for each band, as shown in Table 

HI 

4.3.2. Star clusters 

In order to quantify completeness limits we added artifi- 
cial clusters of different ages and masses to our images. 
We created artificial clusters using a stochastic approach. 
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Table 3. Number of clusters detected per field, per galaxy. Clusters with measured sizes (2^^ column), total of objects 
with three band photometry (3^^^ column), accepted (4^^ and 5^^ columns), suspected (6^^ and 7^^ columns), and total 
rejected (8^^ column). Subscripts 3B and 4B represent three and four band photometry. Shaded areas indicate the total 
per galaxy. 



Galaxy _F 


Ishape 
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Fig. 4. Two-color diagrams for the clusters with four band 
photometry. The red line represents the theoretical path a 
cluster will follow using GALEV models, assuming LMC 
metallicity and no extinction. Accepted^ Suspected clus- 
ters are presented. 

Assuming a Kroupa IMF (|Kroupall2002[ ) in the mass range 
0.01 to 100 Mo, a total cluster mass of M = [10^, 10"^, 10^] 
M0, and a cluster age range between r = [10^,10^*^] yr 
(with 0.5 dex steps), we randomly sampled stars from the 
IMF until the total mass of the stars reached the total 
mass assumed for the cluster. Posi tions were assigned by 
randomly sampling a King profile (|Kingl [19621 ) . There is a 



Table 4. 50 % completeness limits for field stars. 



Galaxy 


Field 


F435W(B) 


F555W(V) 


F814W(I) 


NGC 5236 


1 


26.12 


26.10 


25.10 


NGC 5236 


2 


26.48 


26.35 


25.48 


NGC 7793 


1 


26.60 


26.52 


25.25 


NGC 7793 


2 


26.91 


26.78 


26.03 


NGC 1313 


1 


26.64 


26.66 


26.15 


NGC 1313 


2 


26.65 


26.55 


26.14 


NGC 1313 


3 


26.79 


26.78 


26.41 


NGC 45 


1 


26.77 


26.74 


26.28 


NGC 45 


2 


26.67 


26.60 


26.08 



possible pitfall regarding the mass of the last star sampled, 
because it could overcome the total input (assumed) mass. 
We kept the last star, even if the total mass is higher than 
assumed. This problem affects low-mass clusters more than 
high-mass clusters. With the ages and masses for the stars, 
we then interpolated in isochrones (of LMC-like metallicity) 
from the Padova group ([Marigo et al.ll2QQ8[ ) and assigned 
magnitudes to each star. For all the artificial clusters, an 
FWHM = 2.7 pixels was assumed (corresponding to a 
i?eff ^ 4 pc). Figure [5] shows stamps of artificial clusters of 
different ages and masses, using an average (B, V, and I) 
image of the galaxy NGC 7793 as an example. 

For each combination of age and mass, a total of 100 
randomly generated clusters were added to the science im- 
ages using a square grid. The artificial imag es with clus- 
ters were created using mksynth in BAOLab (.Larsenlll999[ ) 
and added to the science images using imarith in IRAF. 
Following the same procedure used for the cluster pho- 
tometry, an average BVI image was created for each field. 
SExtractor was then run on this average image, using the 
same parameters as in Sect. 4.2. SExtractor returns a file 
with coordinates, measured FWHM and other information. 
To save computational time, owing to the large amount of 
objects that SExtractor could detect, we removed all the 
original objects (science objects detected previously) from 
the list and kept the ones that are not in the original image. 
This new coordinate file was passed to ishape in BAOLab to 
compute PSF-corrected sizes. We used the coordinate file 
to run photometry, again with the same procedures and 
parameters as for the science photometry. Having B, V, I 
photometry done, size cuts were applied and the output 
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Log(M)=3;Log(T)=7 Log(M)=3;Log(T)=8 Log(M)=3;Log(T)=9 




Log(M)=4;Log(T)=7 Log(IVl)=4;Log(T)=8 Log(M)=4;Log(T)=9 




Log(M)=5;Log(T)=7 Log(M)=5;Log(T)=8 Log(M)=5;Log(T)=9 

Fig. 5. Stochastic clusters created to estimate the com- 
pleteness in NGC 7793 (average B, V, and I image, first 
field). From top to bottom, rows represent masses of 
log(M)[M0] = [3,4,5]. From left to right, columns repre- 
sent ages of log(r)[yr] = [7,8,9]. Each images has a size of 
100 X 100 pixels. 



file was matched with the input coordinate file to evalu- 
ate how many of the added artificial objects were recovered 
successfully. 

Figure [6] shows the output (my) average magnitude of 
the recovered clusters versus the fraction recovered for the 
three masses and six ages assumed in each galaxy. For the 
masses log(M)[M0] = [3,4,5] as an illustrative example. Fig. 
[6] shows the total number of objects detected applying three 
different size criteria: (1.) FWHMishape > 0, i.e. no size 
cut used to select the clusters; (2. ) FWHMishape ^ 0.2, 
same size cut used by iMora et al.1 (|2007l |^09); and (3.) 
FWHMishape ^ 0.7, the size cut used in this paper. Each 
line is for a given cluster mass, while each symbol be- 
longing to a line represents the time steps assumed, i.e., 
r = [10^,10^*^] yr (with 0.5 dex step). From this test we 
conclude that 

1. High- mass clusters (log(M)[M0]=5), at any age, are eas- 
ily recognized by our procedures, regardless of the size 
criteria used. 

2. For decreasing cluster mass, clusters of old ages drop out 
of the sample, and for log(r)[yr]>8.5 our completeness 
is less than 50% for masses below log(M)[M0]=4. 

3. The completeness depends on the mass, on the size used 
to classify an object as extended or not extended, and 
on age. Decreasing the size threshold would increase the 




18 20 22 24 26 18 20 22 24 26 18 20 22 24 26 



average V^^^ 

Fig. 6. Completeness curves for stochastic clusters 
with ages log(r)[yr] = [7,7.5,8,8.5,9,9.5] and masses of 
log(M)[M0] = [3,4,5] for the four galaxies. All clusters have 
an input FWHM=2.7 pixels. Red lines represent the num- 
ber of detections without applying any size criteria (right 
column). Blue lines are the number of detections after ap- 
plying the size criteria of FWHMishape > 0.2 pixels, used 
in Mora et al.l (|2007i . 2009) (middle column). Black lines 
are the number of detections after applying the size cri- 
teria used in this work {FWHMishape ^ 0.7 pixels) (left 
column). The legend is in units of solar masses (M0). The 
symbols over the lines represent an age step (i.e. 0.5 dex) 
starting from the left. 

completeness somewhat for low-mass, young objects, 
but could introduce additional contamination. 
4. Stochasticity affects young star clusters (log(r) [yr] <7.5) 
more dramatically. 

At the very young ages in Fig. [6] the completeness drops 
below 100% at all masses. This can be understood from Fig. 
[71 which shows the histograms for the FWHM of the de- 
tected objects measured with ishape for log(M)[M0] = [3,4,5] 
and log(r)[yr] = [7,8,9], using the results from NGC 7793 as 
an example. The figure shows that at high masses and old 
ages, the recovered sizes are on average similar to the in- 
put values, with some spread. However, at younger ages 
and/or lower masses, the measured sizes are systematically 
less than the input values. In these cases, the light profiles 
can be dominated by a single or a few bright stars, while 
for high masses and/or old ages, the light profiles are much 
smoother and better fit by the assumed analytic profiles. 
At young ages and low masses, this bias in the size mea- 
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Fig. 7. Histograms of measured FWHMs for the stochas- 
tic clusters in NGC 7793 created for the completeness 
analysis. From top to bottom, rows represent masses of 
log(M)[M0] = [3,4,5]. From left to right, columns repre- 
sent ages of log(r)[yr] = [7,8,9]. Legends are in logarithmic 
mass and age units. The red dashed line represents an 
FWHM=0.7 pixels, as assumed in Sect. 4.3.2, while the blue 
dash-dotted line represents the input FWHM=2.7 pixels. 



surements leads to a decrease in the completeness fraction 
as more clusters fall below the size cut. 



4.4. Aperture corrections 

Aperture corrections were estimated separately for star and 
clusters. We applied the same procedures as in Paper I. For 
field stars, our PSF-fitting magnitudes were corrected to a 
nominal aperture radius of O'.'S, following standard proce- 
dures. From t his nominal va l ue to infinity, we applied the 
corrections in lSirianni et al.l (|2QQ5[ ). 

For star clusters, ape rture correct i ons w ere applied fol- 
lowing the equations in iMora et aD (|2QQ9[ ) , which give a 
relation between the FWHM of the objects and the aper- 
ture corrections. The photometric parameters (and data 
set) used in our work are the same as the ones used by 
iMora et aD (|2QQ9[ ), allowing us to use their equations. This 
set of equations apply correction s to a nomina l radiu s of 
1''45. We adopted the values in ISirianni et al.l (|2QQ 5^ to 
correct from there on, although the corrections to infinity 
are minor (~ 97% of the total energy is encircled within 



4.5. How do stochastic effects influence star cluster 
photometry? 

The classical approach used to estimate masses and ages of 
unresolved star clusters in the extragalactic field is based on 
the use of multi-band i ntegrated colors and comparison of 
these with SSP models. lAnders et al.l (|2QQ4f ) show that it is 
necessary to have at least four-band photometry to be able 
to break degeneracies (e.g. age-metallicity). However, stan- 
dard SSP models assume a continuously populated stellar 
initial mass function (SIMF), while clusters consist of a fi- 
nite number of stars. For unresolved clusters, the random 
sampling of the SIMF can strongly affect integrated prop- 
erties such as cluster colors, magnit udes, and parameters 
(ages , masses) derived from them (|Cervino fc Luridianal 
l2QQ6l : lMafz ADellanizll2QQ9l : iPopescu fc Hanson 2QlQa.g). A 
promising attempt to take the stochastic color fluctuations 
into acc ount w hen deriving ages and masses has been made 
by iFouesneau fc: LanconI ((2010j), based on a Bayesian ap- 
proach. 

Here we do not attempt to offer any solution to the 
SIMF sampling problem, but we quantify its effects that 
are related to our study. To that aim, we created clusters by 
randomly sampling the SIMF and assigning magnitudes to 
individual stars in the same bands used for our photometry 
(i.e. U, B, V, and I). In addition to its impact on cluster 
detection and classification, as described in the previous 
section, we also investigated how stochastic SIMF sampling 
affects the two-color diagrams and ages. 

We used the same recipe as in Sect. 4.3.2 to create arti- 
ficial clusters. Assuming a range of ages between 10^"^ Myr 
and 10^-^ yr and total masses of M=[10^, 10^ 10"^, 10^ 10^] 
M0, we created 100 clusters every 0.1 dex in age. The top 
row in Fig. [8] shows the two-color diagrams for each one of 
the total masses, together with an solar- met allicity Padova 
SSP model (Marigo et al. 2008) . The evolution of the colors 
U-B and V-I with time are shown in the second and third 
rows, and the comparison of input (assu med) and ou tput 
(estimated) ages using AnalySED (Anders et aDl2004f ) are 
in the bottom row. The colors indicate the input ages. 

Many features are observed here. (1.) For high masses 
(log(Mass)[M0] = [5,6]), the stochastically sampled clusters 
form narrow sequences in the colo r-color a nd color-age dia- 
grams, in agreement with Cervih o fc Luridi ana (2006). (2.) 
At intermediate masses, i.e., log(Mass)[M0] = [3,4], which 
are typical of the cluster masses observed in extragalactic 
works, the scatter increases strongly. The scatter observed 
in the two-color diagrams for this mass range is similar 
to what is seen in our observed two-color diagrams. (3.) 
For very low masses (log(Mass)[M0] = [2]), the scatter again 
decreases but the model colors now deviate strongly from 
the SSP colors. This is because such clusters have a very 
low probability of hosting a luminous (but rare and short- 
lived) post-main sequence star, while the SSP models as- 
sume that the colors are an average over all stages of stellar 
evolution (see Piskun ov et aTl l2009V (4.) For ages log(r)< 8 
and intermediate masses (log(Mass)[M0] = [3,4]), the color 
distribution actual l y beco mes bimodal, as observed by 
iPopescu fc Hansonl (|2010al lbl) . The blue "peak" in the color 
distribution is due to clusters without red supergiants, 
while the presence of even a single red supergiant shifts 
the colors into the other peak. (5.) The bottom row shows 
that age estimates are completely dominated by stochas- 
tic effects for low-mass clusters (log(Mass)[M0] = [2,3]). For 
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Fig. 8. Stochastic effects on colors and ages of clusters. First row.Two color diagrams for a stochastic sample of clusters 
with different masses and different ages. The dashed line represents Padova 2008 SSP models of solar-like metallicity. 
Second row: U-B color evolution. Third row: V-I color evolution. Fourth row: Density plot to make a comparison between 
input and output ages. The red line represents the 1:1 relation, not a fit to the data. 
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higher masses the ages are more accurately recovered, even 
if a small scatter is observed compared to the 1:1 relation, 
especially at ages of a few tens of Myr where the light is 
strongly dominated by red supergiant stars. 

Based on these results it is clear that photometry and 
the ages estimated from broad band photometry can be 
heavily affected by the stochastic effects introduced by 
SIMF sampling, as also shown by Maiz Apellaniz (2009). 
We conclude that low- mass clusters (M < lO^M©) are very 
strongly affected, reaching age differences up to log(agein)- 
\og{digeout) ^ 2.5 dex, for clusters with log(r)[yr]< 8.5, 
while for high-mass clusters (M > lO^M©), this effect grad- 
ually diminishes, an effect that is clearly visible as gaps in 
Fig. [8] (last row). Another important effect is observed as 
a deviation from the 1:1 red line observed for the low-mass 
clusters, which indicate that the estimated age {AgCout) is 
again wrongly recovered, even at old ages (1 Gyr). One 
should be aware of the risk that low-mass, young clusters 
may erroneously be assigned old ages. When using SSP 
models to convert their luminosities to masses, based on 
such wrong age estimates, such objects might be assigned 
erroneously high masses, thus making it into an observed 
sample (e.g. iPopescu fc Hanson|[2010b[ ). 

It is important to mention that the SSP models 
used do not take the binarit y or r otation of mas- 
sive stars into account (s ee e.g. iMaeder fc Mevnet 2008; 
lEldridge fc Stanwavl l2009f ). Also, different isochrone as- 
sumptions and the techniques used to perform the fit to 
the ages and masses can be affecting the res ults (see e.g 
IScheepmak er et all [2QQ1 Ide Griis et ani2005[) . The tests 
presented in this paper are only intended to address the 
stochastic sampling effects, and we only rely on Padova 
isochrones. A more detailed study must be made to better 
account for other effects (e.g. by binaries, rotation, etc). 




0.0 0.5 1.0 1.5 2.0 2.5 0.0 0.5 1.0 1.5 2.0 2.5 

V-l 

Fig. 9. Fitted Hess diagrams for the sample studied. Red 
lines are the boxes used for the fit and the 50% complete- 
ness, same as Fig. [2j 



5. Results 

In the next sections we estimate the field star formation 
histories (SFHs) for each galaxy and determine the ages 
and masses of the cluster candidates. We then estimate the 
cluster formation efficiencies, F. 

5.1. Field star formation histories 

To estimate the SFHs, we used the synthetic CMD method. 
We implemented this method using an IDL-based program 
that was introduced and tested in Paper I. For a description 
of the program we refer the reader to that paper, but here 
we summarize the basic functio nality. 

The synthetic CMD method (jTosi et al.lll99lh takes ad- 
vantage of the power supplied by the CMDs. The method 
uses a group of isochrones, together with assumptions about 
the SIMF, metallicity, distance, extinction, and binarity, to 
reproduce an observed CMD. Photometric errors and com- 
pleteness functions (both being treated as magnitude de- 
pendent parameters by our program) are also taken into 
account. The program searches for the combination of 
isochrones that best matches the observed CMD, thereby 
estimating the SFH. 

The parameters used to estimate the SFH of the galaxies 
are: a Hess diagram with a resolution of 200 x 200 pixels is 
created, using a Gaussian kernel with a standard deviation 
of 0.02 magnitudes along the color axis. The matching is 



done within the rectangular boxes depicted in Fig. [2l using 
the V-I vs. V color comb ination, Padov a 2008 isochrones 
("Marigo et all l2008l and a Kroupa (2002) IMF in the mass 
range 0.1 to 100 M©. The assumed distance moduli, fore- 
ground extinctions, and metallicities are given in Table [TJ 
while the photometric errors and completeness for each 
galaxy were determined in Sect. 4. The program also in- 
cludes a simplified treatment of binaries, in which binary 
evolution is ignored, but the effect of unresolved binaries 
on the CMD are modeled. To account for binarity we used 
three different assumptions for the binary fraction (/) and 
mass ratio {q): (1.) / = 0.0 and q = 0.0, (2.) / = 0.5 and 
q = [0.1,0.9] (assuming a flat distribution), and (3.) / = 1 
and q = 1. These three assumptions are the same ones as 
used in Paper I. 

Figures [9] shows the best-fit Hess diagrams. Comparing 
with the observed Hess diagrams (Fig. [2j), we see that the 
fits are far from perfect. In particular, all the model Hess 
diagrams show a clear separation between the blue core He 
burning ("blue loop") stars and the main sequence, while 
this is not obvious in most of the observed diagrams. This 
might be partly due to some variation in the internal ex- 
tinction, which has not been included in our modeling. To 
infer the SFHs of our sample, we combined all the fields 
(per galaxy) and passed to our program. The star formation 
rates, normalized to unit area, are shown in Figure [TOl and 
the average values are listed in Table [5] for ages between 10 
and 100 Myrs. In this age range, our data are less affected 
by incompleteness. Previous estimates of the Hsfr done by 
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Fig. 10. Star formation rate densities for the data set. Each 
line represent an assumption for the binarity (see text for 
details). 



iLarsen fc Richtlerl (!2000) and ' Chandar et"all (|2QlQb[ ) are 
included in tabled We see that NGC 5236 and NGC 1313 
have higher Ssfr values than NGC 7793 and NGC 45, in 
agreement with the previous estimates. Within the uncer- 
tainties (see Paper I), we do not see any significant trends in 
the SFRs between 10^ and 10^ years. While our estimate 
of SlsFR agrees very well with the others for NGC 5236, 
there are significant differences for some of the other galax- 
ies, most notably for NGC 4395. It should be kept in mind 
that the Ssfr values derived here are for our specific ACS 
fields, while the others are averages over whole galaxies over 
a rather large outer diameter. It is therefore not very sur- 
prising that our new estimates tend to be higher. 

5.2. Cluster ages and masses 

To determine the ages and masses o f the clu sters we 
used the program An alySED ( Anders etap |2QQ4). Using 
GALEV SSP models (|Schulz et a l. 2002), AnalySED com- 
pares the observed spectral energy distributions with a li- 
brary of models to find the best fit. We used GALEV mod- 
els based on a KroupalMF (Kroupa 2002) in the mass range 
0.1 to 100 M0, Padova isochrones (Girard i et al.ll20Q2[ ), and 
different metallicities, depending on the galaxy. 

Based on our analysis of uncertainties due to stochas- 
ticity, we applied a mass criterion to the cluster samples 
in addition to the three selection criteria defined in Sect. 
4. We require clusters in our sample to be more massive 



Fig. 11. Age- mass distributions for the cluster systems. 
Blue dashed lines represent the magnitude cut at rriy = 23. 
Red boxes will be used to create the age and mass distribu- 
tions. Left column present the distribution for the Accepted 
sample, while the right column present the suspected sam- 
ple of clusters (see Sect. 4.2). Red dash-dotted line in the 
right column denotes the mass 10^ Mq for comparison. 



than 1000 M©; however, we remind the reader that ages 
are only reliable for masses greater than 10^ (for ages 
T <l Gyr). The magnitude limit. My ~ —5, is generally 
below our 50% detection limit based on the stochasticity 
test. 

Figure [TTI shows the age- mass diagrams for clusters that 
satisfy the four criteria for the accepted and suspected sam- 
ples separated. Overall, we observe that the suspected clus- 
ters are below 10^ Mq (right column in the figure). We see 
that the number of clusters in NGC 45 is small, compared 
with the other three galaxies. To first order, this appears 
to be consistent with the overall low star formation rate 
derived for this galaxy. Clusters in all four galaxies display 
a range in age and mass, but most are younger than 1 Gyr 
and have masses below 10^ M©, with few exceptions. We 
cannot exclude, however, that the sample contains some 
older clusters that have been assigned too young ages. In 
particular, relatively metal-poor old globulars would not be 
fit well by the models used here. Our sample includes two of 
the eight spectr oscopically c onfirm ed old globular clusters 
in NGC 45 from Mor a et al.l (|2008[ ). From these we find ages 
of 0.04 and 0.8 Gyr, while Mora et al. find spectroscopic av- 
erage ages of 4.5 and 6.5 Gyr but consistent with ages as 
old as 10 Gyr. This confirms the suspicion that some of the 
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Galaxy 


SFRt 

[Moyr-i] 


Areas^ 
[Kpc^] 


[Moyr-^Kpc-^] 


si 

[Moyr-^Kpc-^] 


[Moyr-^Kpc-^] 


NGC5236 


0.39 


28.71 


13.43 X 10"^ 


13.8 X 10"^ 


16.8 X 10"^ 


NGC7793 


0.15 


23.05 


6.43 X 10"^ 


2.12 X 10"^ 


2.8 X 10"^ 


NGC1313 


0.68 


60 


11.26 X 10"^ 


4.04 X 10"^ 




NGC45 


0.05 


48.99 


1.01 X 10"^ 


0.23 X 10"^ 




NGC4395^ 


0.17 


36.48 


4.66 X 10-^ 


0.25 X 10"^ 





estimated in this paper; ^Larsen k, Richtle^ (|200QD : " iChandar et al.l (|2010br ): ^ values from Paper I. 



clusters in our sample might be older, metal-poor globulars 
with misassigned ages. 

Our catalogs of the cluster candidates are available on- 
line. Table [10] shows a few example lines to illustrate the 
format and the information contained in the catalogs. 

5.3. Cluster disruption and formation efficiencies 

Although we have derived masses and ages for our cluster 
candidates, some additional steps are necessary before we 
can use this information to der ive cluster for mation rates. 
In Paper I we used a Schechter (|Schechterlll976 ) mass func- 
tion (M^ = 2 X 10^ M0) to extrapolate below the (age- 
dependent) mass Hmit and in this way we estimated the 
total mass in clusters with ages between 10^ and 10^ years 
for NGC 4395. This approach ignores any effects of disrup- 
tion but is still useful for relative comparisons. We therefore 
first apply the same approach to the four galaxies in this 
paper. Completeness limits were estimated by plotting the 
luminosity functions and identifying the point where they 
start to deviate significantly from a smooth power law. This 
occurs at the following absolute magnitudes: My = —6.2 
for NGC 5236, My = -5.7 for NGC 7793, My = -6.8 for 
NGC 1313, and My = -5.9 for NGC 45. After estimating 
the mass in clusters with 10^ < r/yr < 10^ down to a limit 
of 10 M0 and dividing by the age interval (see Paper I for 
details), the resulting CFRs were normalized to the area of 
the full ACS fields (a factor of ~ 2.27 more) for comparison 
with the field star formation rates. The resulting CFRs are 
listed in the second column of Table[8] (CFRpi). 

Especially for NGC 45, the CFRs derived in this way 
are highly uncertain owing to the small number of clusters 
that have four-band photometry. Better statistics can be 
obtained by only using the three-band photometry in the 
ACS frames, but at the cost of having no age information for 
individual clusters. However, CFRs may still be estimated 
by comparing the observed luminosity functions (LFs) with 
scaled model LFs (Gieles 2010). If the CFR is assumed 
constant and assumptions made about the initial cluster 
mass function (^) and disrupti on pa rameters, the LF can 



on pa 

be modeled as follows (Eq. 7 in lLarsen.2QQ9.) : 



dN_ _ r 



^,[M,(L,r)] X 



dMc 

x/surv(T) dr , 



X Tc(r) X CFR 



(1) 



where ^i[M^(L, r)] is the initial mass function; Tc(r) is the 
mass-to-light ratio, which is only dependent on time in or- 
der to be able to compute it from classical SSP models; CFR 
is assumed constant over time; and /surv(T) is the number 
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Fig. 12. Age distributions for clusters in the galaxies 
NGC 5236, NGC 7793, and NGC 1313 in the mass range 
10^ to 10^ M0. Black lines represent the distributions for 
galaxies. Red lines represent the best fit. Upper panel are 
the ADs only using the clusters classified as accepted. The 
lower panel are the ADs for the clusters classified as ac- 
cepted plus suspected. 



of clusters that survive after applying MID; Mi and Mc are 
the initial and current masses of a cluster with luminosity 
L and age r, and these are related through the assumed 
secular disruption law (Lamers et al. 2005). The descrip- 
tion of MID adopted here assumes that a constant fraction 
of the cluster population is removed per logarithmic age 
bin, as opposed to a constant fraction of the mass of each 
individual cluster. If is a uniform power law, this makes 
no difference. However, if the MID stem from gradual mass 
loss from individual clusters, any features in (such as 
the cut-off mass for the Schechter function) will be shifted 
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Table 6. Slopes for the age distributions for the accepted Table 7. Slopes for the age distributions for the accepted 



and the accepted+suspected samples. 



Mass [Mo] 


NGC 5236 NGC 7793 


NGC 1313 


10^ - 10^ 


Accepted 
-0.23±0.1 -0.72zb0.27 


-0.63±0.12 


10^ - 10^ 


Accepted+Suspected 
-0.26±0.08 -0.42±0.19 


-0.57±0.11 



and the accepted+suspected samples based on the maxi- 
mum likelihood fit. 



Mass [Mo] 


NGC 5236 NGC 7793 


NGC 1313 


10^ - 10^ 


Accepted 
-0.27±0.11 -0.48±0.32 


-0.63±0.10 


10^ - 10^ 


Accepted+Suspected 
-0.31=b0.09 -0.30=b0.21 


-0.60=b0.09 



downwards with time, whereas only the normalization of 
will change with time for constant number loss. 
In order to apply Eq. ([T]), some constraints on clus- 
ter disruption are necessary. Models and empirical con- 
straints on cluster disruption have been discussed in re- 
cent y ears by different authors (e.g 



2003[ iLamers et al l 120051 : IWh 



Boutlou. 



2009[ IZhang fc Falll [l999l : iFalll l2004 among others), and 



it more et al.l 



cos fc 
20071 : 



jamers 



Larsen 



Fall et al.l (|2009f )"for different types of galaxies such as the 



LMC, SMC, Milky Way, M83, or Antennae. Given that it 
is currently uncertain to what extent MID or MDD domi- 
nates the cluster disruption, we carried out our analysis for 
both scenarios. 

Figure [12] shows the age distributions (ADs) for clus- 
ters with masses between 10^ and 10^ M© for the galaxies 
NGC 1313, NGC 5236, and NGC 7793. NGC 45 has too 
few clusters to derive meaningful ADs. We show fits for 
both the Accepted and Accepted+Suspected samples. The 
slopes of the age distributions, obtained by carrying out 
fits of the form \og{dN/dt) = a x log(r) -\- b to the data 
in Fig. [121 given in Table [6] There are no large dif- 
ferences between the slopes derived for the Accepted and 
Accepted+Suspected sample. Figure [TTI shows that most of 
the clusters in the suspected sample have masses below our 
limit of log(M) [M©] = 4.0, explaining the similarity of the 
age distributions above this limit. 

As a consistency check for the slope of the age distribu- 
tions, we performed a maximum likelihood fit to the data, 
assuming a power-law relation and using the power-law in- 
dex as a free parameter. Using the accepted and accepted 
plus suspected sample of clusters, we estimated the slope of 
the age distributions using the same age and mass ranges 
as shown in Fig. [12] (i.e. ages between 4 Myrs up to 1 Gyr 
and masses between 10^ and 10^ M©). The results obtained 
are presented in Table [7] The derived slopes agree very well 
with those in table [6] within the errors. 

Using clusters with ages between 10^*^ < t < 10^ yr 
and a mass 10"^ < M < 10^ Mq we checked for (possible) 
variations over the slope of the age distributions. The best 
fit for these slopes in the new age interval are —0.17 + 0.39, 
—0.38 + 2.34, and 0.78 + 0.43 using the accepted sample, 
and -0.39 + 0.31, 0.05 + 0.91, and 0.59 + 0.91 using the 
accepted plus suspected sample, for the galaxies NGC 5236, 
NGC 7793, and NGC 1313 respectively. The new slopes 
are to be flatter than the values for the whole age range, 
showing even positive values; however, in most of the cases, 
the error is larger than the estimation itself. The possibility 
of having shallower slopes indicate that there is a possible 
curvature of the age distribution of star cluster systems, 
which is not consistent with a (simple) power law. 

The slopes found here are, however, somewhat shallower 
than the value of a = -0.9 + 0.2 found for NGC 5236 by 
IChandar etliD (|2010cD . If interpreted within the MID sce- 



nario, the slopes derived here correspond to MID disruption 
rates of 41%, 81%, and 77% per decade in age if we use the 
accepted sample, and 45%, 62%, and 73% if we use the 
sample of accepted plus suspected objects, for the galax- 
ies NGC 5236, NGC 7793, and NGC 1313, respectively. A 
weighted average of the slopes for the Accepted samples 
leads to a mean slope of (a) = —0.42 + 0.07 and an MID 
disruption rate of (62 + 6)% per decade in age. We use 
this mean weighted value for all the galaxies in our sample, 
including NGC 45 and NGC 4395 where the numbers of 
detected clusters are too low to allow us to estimate the 
slopes of the age distributions independently. 

It is wor th comparing our age distributions with those of 
iM ora et aD ([2009). Their slopes were considerably steeper 
than those found here, but this is partly because Mora et al. 
worked with magnitude-limited samples, while our age dis- 
tributions are for mass-limited samples. When taking this 
into account. Mora et al. found that their data were consis- 
tent with MID disruption rates of 75%-85% per dex, still 
somewhat higher than the values found here. We note that 
the stochastic effects, combined with our relatively conser- 
vative size cuts, likely cause us to underestimate the num- 
ber of objects in the youngest bins. This might account 
for some of the flattening of the age distributions in the 
youngest bins that is also seen in Fig. [12] for NGC 1313 
and NGC 5236. For this reason our estimates of the slopes 
and disruption rates might also be considered lower limits. 
The less strict size cut used by Mora et al. (FWHM=0.2 
pixels instead of our 0.7 pixels) would cause them to detect 
more compact, young objects, but also possibly ones with 
increased contamination. Finally, since visual selection was 
part of the sample selection in both our work and the one 
of Mora et al., this will also lead to differences in the final 
samples. 

A key difference between the MID and MDD models is 
that, while the former predicts no change in the shape of 
the mass distribution with age, the latter predicts a flat- 
tening at low masses. In order to test whether we can dis- 
tinguish between the two scenarios, model mass distribu- 
tions (MD, number of objects, per mass bin, per linear 
age bin) were computed for different age intervals using 
a relation analogous to Eq. [TJ but integrating over r for 
fixed M rather than fixed L. These model MDs were then 
compared with the observed MDs in the same age inter- 
vals. For the MDD scenario we assumed 7 = 0.62 and a 
disruption time of ^4 = 1 x 10^ yr (jLamers et al.l 120051 ). 
No infant mortality was included in the MDD models (i.e. 
/surv = 1 at all ages). The MID models used the IMR ob- 
tained above (62%), which is active in the age range 5 Myrs 
to 1 Gyr, and was set to infinity, making the models in- 
dependent of mass. We compared the model MDs with our 
observations, dividing our cluster samples into the follow- 
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Fig. 13. Mass-independent models. Left column: MDs ob- 
served (black symbols) and predicted models (red lines). 
Right column: Luminosity function observed (black hori- 
zontal lines), predicted models (red horizontal lines), and 
theoretical model (dash-dotted line). Vertical straight lines 
represent the limits used for the fit of the LFs. 



Fig. 14. Mass-dependent models. Left column: MDs ob- 
served (black symbols) and predicted models (red lines). 
Right column: Luminosity function observed (black hori- 
zontal lines), predicted models (red horizontal lines), and 
theoretical model (dash-dotted line). Vertical straight lines 
represent the limits used for the fit of the LFs. 



ing age and mass bins: log(M)[M0] = [3.O,3.5,4.O,4.5,5.O] and 
log(r)[yr] = [6.6,7.5,8.0,8.5,9.0] (shown in Fig.^. 

In Figs. [13] and [14] we compare the observed and pre- 
dicted MD and LF for MID and MDD models. The best 
fits were defined by scaling the model LFs to match the 
observed ones, where the observed LFs were created using 
variable binning with ten objects per bin for N GC 1313 
and NGC 5236 f following iMaiz Apellanizl (|2QQ9[ )) and five 
objects per bin for NGC 45, NGC 4395, and NGC 7793. 
Errors were assumed to be Poissonian. Also shown is the 
reduced foi" each galaxy. Based on these figures, it is dif- 
ficult to distinguish between the MID and MDD scenarios, 
mostly because of the limited dynamical range of the data 
and the poor statistics. 

Our final estimates of the CFRs come from the LFs 
for clusters with 3-band photometry in all five galaxies in 
our sample, using model LFs computed for both MDD and 
MID disruption. The observed LFs were created in the same 
way as described above, using variable binning. We used 
the magnitude ranges up to the brightest cluster in the 
sample and the lower limit was set at the completeness 
limits estimated above (see vertical lines in Fig. 15). The 
CFRs are listed in Table 7 for the 4- and 3 band photometry. 
The CFRs for the four-band photometry have been scaled 
to the size of the ACS fields. The observed LFs and the best- 
fitting models are shown in Fig. [151 where MID and MDD 



models have been tested for the five galaxies in our sample. 
The CFRs inferred from the fits in Fig 15 only include 
the accepted clusters. If we include the "suspected" objects 
then the CFRs increase by 52%, 62%, 30%, 4%, and 25% for 
NGC 5236, NGC 7793, NGC 1313, NGC 45 and NGC4395, 
respectively, for the MDD scenario. Similar changes occur 
for MID. Fur thermore, if we e xclude the known ancient GCs 
in NGC 45 ( Mora et al.ll2QQ8[ ) from the sample, the F value 
for this galaxy decreases by ^ 20%. Taken together, this 
then makes NGC 45 less of an outlier in Fig. [16] and in 
general shifts the data points upwards. 

Having estimated the SFRs and CFRs, we are now able 
to calculate F. Table [9] lists the F values corresponding to 
the different estimates of the CFR. As an additional con- 
sistency check, we also include a direct comparison of the 
fraction of light coming from clusters and from field stars, 
using the magnitude range between my = 23 and my = 18. 
This is not a direct measure of F, but is still useful for check- 
ing any trends. The F values derived from the LF fitting 
are generally a few percent, with the values derived from 
MDD models about 48% smaller (or less) than those de- 
rived for the MID models. We computed the mean value 
for each model and the respective standard error of the 
mean (columns 5th and 8th in Table [9]). 
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Table 8. Estimates of the cluster formation rates. Subscript PI refers to estimations made following Paper I. Disruption 
models are labeled as MDD and MID. Number of bands used for the estimation are labeled as 3B and 4B. 



Galaxy 


CFRpi 

[Moyr-i] 


K^r iXMDD 

[Moyr-i] 


^MDD 

[Moyr-i] 


^MID 

[Moyr-i] 


^MID 

[Moyr-i] 


NGC5236 


37.7x10"^ 


23.0x10"^ 


19.9x10"^ 


41.0x10"^ 


36.6x10-^ 


NGC7793 


14.5x10"^ 


4.0x10"^ 


3.4x10-^ 


6.8x10-^ 


6.6x10-^ 


NGC1313 


60.7x10"^ 


23.0x10"^ 


19.7x10"^ 


44.0x10"^ 


38.6x10"^ 


NGC45 


8.6x10"^ 


2.5x10"^ 


2.7x10"^ 


4.4x10"^ 


3.6x10-^ 


NGC4395 


4.5x10"^ 


2.4x10"^ 


1.0x10"^ 


4.0x10"^ 


1.6x10"^ 



Table 9. Estimations of F and comparison with iGoddard et aD (|2Q1Q| ) results. 



Galaxy 




p3i^ 
J- MDD 


J- MDD 


^MDD ± CFMDD 


p3i^ 
-L MID 


-L MID 


^MID ± CFMID 


rt 


TL{uy 


Lcl/Lfs 




[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 


[%] 




[%] 


NGC5236 


9.8 


5.9 


5.2 


5.6±0.6 


10.5 


9.4 


10.0=b0.9 


10.3 


2.36±0.31 


9.0 


NGC7793 


9.8 


2.6 


2.3 


2.5±0.3 


4.5 


4.5 


4.5=b0.1 


8.6 


1.15=b0.32 


3.4 


NGC1313 


9.0 


3.3 


2.9 


3.2±0.2 


6.4 


5.7 


6.1±0.6 


9.9 


1.49zb0.44 


11.7 


NGC45 


17.3 


5.0 


5.4 


5.2±0.3 


8.8 


7.2 


8.0=bl.l 


2.2 


0.24zb0.17 


6.1 


NGC4395 


2.6 


1.4 


0.6 


1.0±0.6 


2.3 


0.9 


1.6±0.1 


8 


0.07±0.05 


1.7 



iLarsen Richtle^ ([2000') . Super and sub- indices for V are the same as Table O Columns 5 and 8 are the mean and standard 
error for each model. Column 9 values using Goddard et al. with our T^sfr- Column 11 is the fraction between the light from 
clusters and the light from field stars in the magnitude range 18 to 23 using My- 



6. Discussion 

In Fig. [T6lwe compare our F vs. Ssfr measurements with 
the data in lGoddard et al.] (|2010f ) (rhombs). We cannot con- 
firm whether a correlation is present within the range of 
^SFR values probed by our data, but overall our F values 
are similar to those found by Goddard et al. in this range 
of HsFR or slightly low er. However, f ollowing the methods 
of Goddard et al., Adan iQ~et al.l (|201 1) estimated values for 
F in the two blue compact galaxies ESO 185 and Haro 11. 
The results from Adamo et al. are in good agreement with 
the power law proposed by iGoddard et al.l (|2010h . 

According to iLarsen fc Richtlerl (|2000 ^ the five galax- 
ies span a significant range in specific [/-band luminosity. 
From Tableland Fig.[T71 we see that galaxies with high F 
values generally tend to also have high Tl{U) values. One 
exception to this is NGC 45, which has a rather high F 
for its Tl{U). The Tl{U) value for this galaxy is, however, 
based on only two clusters, hence subject to very large un- 
certainty. Nevertheless, the high F value for NGC 45 is also 
somewhat puzzling given that it has the lowest Ssfr- This 
may suggest that there is not a simple relation between F 
and SsFR- In this context, it is also interesting that this 
galaxy has a large number of ancient GCs for its lumi- 
nosity, yielding an unusually high glo bular cluster spe cific 
frequency for a late- type (Sd) galaxy ('Mora et al.l l2QQ9l) . 

Our measurements of F values in the range ^2-10% 
are consistent with other recent estimates of the fraction 
of stars forming in bound clusters. It should be kept in 
mind that this number is not necessarily an indicator of 
"clustered" vs. "isolated" star formation, since some stars 
may form in embedded clusters that dissolve or expand on 
short enough time scales to drop out of our sample. 

There is a correlation between Tl(U) and F, as shown 
in Fig. [T71 We estimated values for F using the ac- 
cepted+suspected samples for the three-band photometry 
and for MDD and MID models as presented in the figure. 



The galaxy NGC 45 deviates from the apparent relation. 
Two things must be noted. (1) The inclusion of the sus- 
pected sample does not change the trend observed dramat- 
ically, although the increases in the CFRs is reflected in 
the new F estimations, and (2) CFR estimates based on 
different disruption models follow the same trend. 

Estimates of the actual "infant mortality rate" are 
hard to make unless the embedded phase is probed di- 
rectly, some thing which is d ifficul t in external galaxies. For 
NGC 1313, iPellerin et all (|2QQ7D find that the IMR is a 
very efficient process for the dissolution of star clusters 
in this galaxy {IMR = 90%) based on UV fluxes in and 
out of clusters. Our estimate of a F value of 3%-5% for 
NGC 1313 indicates that ^ 95% of star formation in this 
galaxy happens outside clusters that are detected in our 
sample, in reasonable agreement with the Pellerin et al. es- 
timate. However, it is also clear that cluster dissolution is 
a continuous process, and systems that are probed at older 
ages are generally expected to show a lower fraction of stars 
in clusters. 

A proper account of dissolution effects could, in prin- 
ciple, be used to correct measurements of F at different 
ages to a common reference (say, 10 Myr), but current un- 
certainties in the disruption process mak es this difficult 
to apply in practice. As a case in point, IChandar et al.l 
(2010a) estimate a mass-independent disruption rate of 
80%-90% per decade in age for NGC 5236, a value that 
is significantly higher than our estimate ~ 40% (~ 62% 
is our estimated weighted average). These differences un- 
derscore that the definition of cluster samples (especially 
in star-forming galaxies) is subject to strong selection ef- 
fects, many of which are age dependent and thus likely to 
affect the age distributions. One example is the size cuts, 
which can easily cause a bias against young objects where 
stochastic SIMF sampling leads to underestimated sizes 
even for masses of ~ lO^M©. Furthermore, there may also 
be a physical relation between cluster size and age (e.g. 
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Fig. 17. Relation between Tl(U) and F using the 
models MDD and MID. The values for T using the 
Accepted+Suspected sample are also included. 



When using young clusters (r less than 100 Myr) we ob- 
served that the slope of the age distribution gets shallower, 
indicating a possible curvature or deviation from a power- 
law relation and showing values for the age distributions 
different from previous estimations. However, these results 
are based on a sample that is strongly affected by very few 
of clusters. 



Fig. 15. LFs for the five galaxies in our sample. We used 
MDD and MID theoretical models in left and right columns, 
respectively. Vertical lines are the limits of the fit. Red hor- 
izontal lines represent the binned theoretical values, while 
black horizontal lines represent observations. Errors are 
Poissoninan. Dashed-dotted line is the theoretical LF not 
binned. 
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Fig. 16. Relation between F and T^sfr using average values 
of F for the models MDD and MID . Rho mbs symbols and 
line represents the iGoddard et al.l (|2Q10[ ) data and black 
star symbols represent our set of galaxies. 



Elson et al.l 119891 : iBarmbv et all 120091 : iMackev fc Gilmord 
2003al lbl). compounding this problem. 



7. Summary and conclusions 

Using HST observations of the galaxies NGC 45, NGC 1313, 
NGC 4395, NGC 5236, and NGC 7793, we studied their 
populations of star clusters and field stars separately with 
the aim of constraining the quantity F, i.e. the ratio of 
stars forming in bound clusters and the "field". We have 
been following the basic approach described in Paper I, 
i.e. comparing synthetic and observed color-magnitude di- 
agrams (for the field stars) and SED model fitting (for the 
star clusters), to get the formation histories of stars and 
clusters. 

We tested how stochastic effects induced by the SIMF 
influence photometry and the estimation of ages and how 
the completeness limits are affected. We conclude that mas- 
sive clusters (log(Mass)[M0] > 5]) are easily detected (with 
the parameters used in this work) at any age, while the 
detection of clusters with masses below - 10"^ Mo can be 
strongly affected by stochasticity. Our tests thus show that 
completeness functions do not just depend on magnitude, 
but also on age and size. It would be desirable to find better 
classification methods than a simple size cut to determine 
what is, and what is not, a cluster. 

We estimated star formation histories and found that 
NGC 5236 and NGC 1313 have the highest star formation 
rates, while NGC 7793, NGC 4395, and NGC 45 have lower 
SFRs. Within the uncertainties, we do not see significant 
variations within the past 100 Myr. 

Comparing observed and modeled mass- and and lu- 
minosity distributions for the cluster populations in differ- 
ent galaxies, we find that we cannot distinguish between 
different disruption models (mass dependent vs. mass in- 
dependent). We compared model luminosity functions for 
each disruption scenario with observed LFs and derived 
CFRs for the cluster systems. From our measurements of 
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the CFRs and SFRs we derived the ratio of the two, F, as 
an indication of the formation efficiency of clusters that re- 
main identifiable until at least 10^ years. We find F values 
in the range ~2-10%, with no clear correlation with Ssfr 
within the (limited) range probed by our data. However, 
our measurements ar e roughly consistent with those of 
iGoddard et al.l (|2Q1Q[ ), who find a relation between Ssfr 
and F for a sample of galaxies spanning a wider range in 
^SFR (but more heterogeneous data). 

A general difficulty in this type of work is to identify 
a reliable sample of bona-fide clusters. Comparison with 
previous work suggests that the cluster samples, although 
covering the same galaxies, are significantly affected by the 
criteria used to classify clusters over the images. This re- 
sults in different estimates of cluster system parameters, 
such as those related to the disruption law. Accurate esti- 
mates of these parameters are also hampered by the rela- 
tively poor statistics that result from having only patchy 
coverage of large, nearby galaxies in typical HST imaging 
programs. 
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Table 10. Online material. The rows presented in this table illustrate what the online material will look like. Column 1st: Name of the galaxy, field observed and 
cluster number. Column 2nd: ID classification .Columns 3rd and 4th: X and Y coordinates of the clusters over the images. Columns 5th and 6th: Right ascension 
and declination (J2000). Columns 7th- 14th: U, B, V, and I magnitudes and their respective errors. Column 15th: Age estimated for LMC-like metallicity. Column 
16th: Mass estimated for LMC-like metallicity. Column 17th and 18th: Sizes measured with SExtractor and Ishape, respectively. Column 19th: Flag for accepted 
and suspected objects. 



Galaxy _Field_# 


ID 


X 


Y 


RA 


DEC 


U 


Ue 


B 


Be 


V 


Ve 


I 


le 






[FIX] 


[FIX] 


[J2000] 


[J2000] 


[Mag] 


[Mag] 


[Mag] 


[Mag] 


[Mag] 


[Mag] 


[Mag] 


[Mag] 


NGC5236_1_1 


5214 


245 


2539. 


13:36:58.81 


-29:51:13.94 


21.08 


0.10 


20.89 


0.01 


20.86 


0.01 


20.57 


0.03 


NGC5236_1_2 


5238 


340 


2549. 


13:36:58.99 


-29:51:09.77 


22.90 


0.21 


22.27 


0.02 


21.66 


0.02 


20.64 


0.02 


NGC5236_1_3 


4360 


596 


2225. 


13:36:58.23 


-29:50:51.62 


22.55 


0.15 


22.32 


0.02 


21.90 


0.02 


21.09 


0.02 


NGC5236_1_4 


5454 


614 


2609. 


13:36:59.62 


-29:50:58.36 


20.73 


0.04 


20.64 


0.01 


20.54 


0.01 


20.14 


0.01 


NGC5236_1_5 


3919 


621 


1928. 


13:36:57.22 


-29:50:44.62 


21.46 


0.07 


21.79 


0.02 


21.76 


0.02 


21.24 


0.03 


NGC5236_2_1 


6116 


574. 


2690. 


13:37:05.64 


-29:56:51.97 


22.07 


0.08 


22.15 


0.01 


22.04 


0.01 


21.61 


0.02 


NGC5236_2_2 


5845 


655. 


2423. 


13:37:04.84 


-29:56:42.65 


22.05 


0.13 


21.65 


0.02 


21.40 


0.01 


20.71 


0.02 


NGC5236_2_3 


5690 


719. 


2180. 


13:37:04.10 


-29:56:34.62 


22.29 


0.12 


22.57 


0.02 


22.28 


0.02 


21.62 


0.02 


NGC5236_2_4 


5859 


746. 


2461. 


13:37:05.12 


-29:56:39.33 


21.72 


0.08 


21.35 


0.01 


21.17 


0.01 


20.71 


0.01 


NGC5236_2_5 


5773 


821. 


2297. 


13:37:04.67 


-29:56:32.47 


21.87 


0.08 


21.97 


0.01 


21.74 


0.01 


21.23 


0.02 



Table 11. Continuation Table [TOl 



Galaxy _Field_# 


Log(r) 


Log(M) 


FWHM^^;, 


FWHMz,;,, 


Flag 




[Yr] 


[Mo] 


[Fix] 


[Fix] 




NGC5236_1_1 


8.00 


3.86 


7.35 


2.91 


acpt 


NGC5236_1_2 


8.87 


4.43 


5.84 


0.74 


acpt 


NGC5236_1_3 


8.92 


4.08 


3.81 


1.30 


acpt 


NGC5236_1_4 


8.39 


4.20 


4.80 


1.93 


acpt 


NGC5236_1_5 


8.05 


3.53 


6.16 


1.63 


acpt 


NGC5236_2_1 


8.00 


3.48 


3.89 


0.95 


acpt 


NGC5236_2_2 


8.78 


4.16 


7.37 


3.67 


acpt 


NGC5236_2_3 


7.81 


3.48 


3.91 


1.35 


acpt 


NGC5236_2_4 


8.70 


4.18 


6.79 


2.79 


acpt 


NGC5236_2_5 


8.39 


3.73 


4.39 


1.84 


acpt 



